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1 INTRODUCTION 

According to the paradigm that the X-ray background 
(XRB) originates mainly from sources at redshifts 1 < 
z < 3 (eg Shanks et al. 1991), it should provide a means 
of measuring the well established solar motion with re- 
spect to the cosmic rest-frame, defined by the Cosmic 
Microwave Background (CMB), towards I — 264°, 6 = 
48° (cf. Lineweaver et al 1996) . An imprint of this reflex 
motion should be a dipole pointing towards the direc- 
tion of the CMB dipole (Compton-Getting effect). The 
available all-sky X-ray maps (HEAO-1 and ROSAT ) 
are composed by X-ray counts not only originating from 
such distant sources but also from local extragalactic 
sources (z < 0.1). These gravitating extragalactic ob- 
jects, which cause our peculiar motion with respect to 
the cosmic rest frame, emit X-rays and therefore their 
dipole should also point towards the CMB dipole. This 
has been demonstrated for the case of AGNs (Miyaji 
& Boldt 1994) and of X-ray clusters (Lahav et al 1989; 
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We estimate the dipole of the diffuse 1.5 keV X-ray background from the 
ROSAT all-sky survey map of Snowden et al (1995). We first subtract the 
diffuse Galactic emission by fitting to the data an exponential scale height, 
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respect to the cosmic rest frame (Compton-Geting effect) we find the large- 
scale dipole of the X-ray emitting extragalactic sources having an amplitude 
2^lss ~ 0.92?xrb, in general agreement with the predictions of Lahav et al 
(1997). We finally estimate that the Virgo cluster is responsible for ~ 20 per 
cent of the total measured XRB dipole amplitude. 
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Plionis & Kolokotronis 1998). Therefore, the dipole pat- 
tern of the XRB results from at least two effects; (a) the 
motion of the observer with respect to the XRB which in 
the context of the cosmic-ray background was discussed 
first by Compton & Getting (1935) and (b) the X-ray 
emission of extragalactic objects the gravitational field 
of which causes the observers motion with respect to the 
cosmic rest frame. The relative contribution of these ef- 
fects have been analytically estimated by Lahav, Piran 
& Treyer (1997) and were found to be of the same order 
of magnitude. 

The dipole anisotropy of the XRB has been 
measured in hard X-rays (2-10 keV) using UHURU 
(Protheroe, Wolfendale & Wdoczyk 1980) and HEAO- 
1 data (Shafer & Fabian 1983, Jahoda 1993). The result- 
ing dipole was found to point towards the general direc- 
tion of the CMB dipole but with a large uncertainty: the 
90 per cent confidence levels quoted by Shafer & Fabian 
(1983) cover 12 per cent of the whole sky. The use of the 
ROSAT all-sky survey can extend these studies to soft 
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energies with higher sensitivity and angular resolution. 
Kneissl et al. (1997) presented a cross-correlation of the 
COBE DMR and the ROSAT all-sky survey maps. How- 
ever, their attempt to measure the extragalactic dipole 
was hampered by the contamination of the Galactic 
emission which is appreciable, especially at low Galac- 
tic latitudes, even in the hard ROSAT band. They con- 
cluded that proper modelling of the Galactic contribu- 
tion is necessary in order to obtain a measurement of the 
dipole. In this letter, we attempt to model the diffuse 
Galactic component using a simple disk model. After 
subtracting our Galaxy model, we estimate the cosmo- 
logical X-ray dipole using the ROSAT all-sky survey 
hard maps (1.5 keV) of Snowden et al. (1995). 



2 THE ROSAT ALL-SKY SURVEY DATA 

We use the ROSAT all-sky survey maps at a mean en- 
ergy of 1.5 keV (PSPC PI channels 91-201). These maps 
are now publicly released and are described in detail in 
Snowden et al. (1995). The maps cover ~ 98 per cent 
of the sky with a resolution of 2 degrees. Point sources 
detected in the all-sky survey are included in the maps. 
The particle background, scattered solar X-ray back- 
ground and other non-cosmic background contamina- 
tion components have been removed (see Snowden et 
al. 1994, Snowden et al. 1995). 

The X-ray emission in the 1.5 keV band is mainly 
extragalactic: Hasinger et al. (1998) have resolved about 
70 per cent of the background at these energies into 
discrete extragalactic sources. However, there is some 
contamination, due to the poor energy resolution of 
the ROSAT PSPC, from a Galactic component. This 
is well-fitted with a Raymond-Smith spectrum, having 
a temperature of ~0.17 keV, and may be associated 
with the Galaxy halo (see Wang & McCray 1993, Gen- 
dreau et al. 1995, Hasinger 1996, Pietz et al. 1998). 
We note that at higher energies (3-60 keV) there is 
evidence for a even harder Galactic component with a 
bremsstrahlung spectrum of 9 keV (Iwan et al. 1982); 
however, this is expected to contribute less than one 
per cent in the ROSAT 1.5 keV band. Moreover, the 1.5 
keV maps of Snowden et al. (1995) show some extended 
features superimposed on the extragalactic and the dif- 
fuse hard Galactic component (see Snowden et al. 1995) 
mainly originating from nearby supernova remnants (eg 
the North Polar Spur, the Cygnus superbubble) . All the 
above local features must be subtracted before deriving 
the extragalactic X-ray dipole. 



2.1 Modelling the Galactic emission 

The derivation of a detailed Galactic emission model re- 
quires observations in many wavebands and is outside 
the scope of this paper; detailed modelling of the Galac- 
tic halo is discussed in Iwan et al. (1982), Nousek et al. 
(1982) and Pietz et al. (1998). Here instead, we attempt 
to make a rough model of the Galactic contamination in 
our energy band. We model the diffuse Galactic compo- 
nents with a finite radius disk with an exponential scale 
height, (eg Iwan et al. 1982) which provides a good de- 
scription of the Galactic component at both soft (0.75 
keV) and hard energies (2-60 keV). This is given by: 



C{l,b)=Cb 



1 + 



Eh 
sin |6| 



-f(l,R d )t a n\b\/h 



(1) 



where f(l, Rd) = cos / + \J R? d — sin 2 I, with C the total 
X-ray intensity, Cb the average extragalactic component, 
in units of 10 -6 cts s _1 arcmin~ 2 , E the fraction of the 
total X-ray emission which is due to the Galaxy, h and 
Rd the disk scale height and disk radius, respectively 
both in units of 10 kpc (the galactocentric distance of 
the Sun). We exclude from the fit the regions of the most 
apparent extended emission features: the bulge and the 
North Polar Spur (ie. -40° < b < 75° and 300° < 30° 
eg Snowden et al. 1997) as well as the Galactic plane 
strip (with \b\ < 20° or 30°). Although it is unknown 
whether there is some small residual bulge emission out- 
side the above excised region, our rough model should 
provide a good first order approximation to the Galac- 
tic halo emission. Furthermore, applying a homogeneous 
mean count we mask the most apparent "local" extra- 
galactic features; a 4° radius region around the Virgo 
cluster (!,!)« 287°, 75°) and a 10° radius region around 
the Magellanic clouds (/, b « 278°, -32°). 

In the minimization procedure we weight each pixel 
by 1/cr, where a = (of + erf,) 1 / 2 , with erf the vari- 
ance of the X-ray ROSAT counts due to the intrinsic 
extragalactic fluctuations and erf, the variance due to 
Poisson count statistics. We estimate 07 excluding from 
the map the North Polar Spur and the |6| < 45° re- 
gions and find 07 ~ 0.27C Starting from different ini- 
tial guesses of the input parameters the minimization 
procedure does not reach a unique minimum, although 
the reduced \ 2 ls around one and the output model 
parameters are closely clustered. This suggests the ex- 
istence of a broad and shallow minimum. We therefore 
run 1000 x 2 minimizations starting from a broad range 
of initial values. These values are centred on those that 
Iwan et al (1982) find for the Galactic component us- 
ing HEAO-1 data, ie., (h) « 7 kpc, (R d ) w 28 kpc and 
10 per cent; while the input extragalactic contri- 
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Figure 1. Distribution of the finite Galaxy disk model parameters. Continuous lines represent the output parameter distribution 
while the hatched histograms the corresponding input ones, used to start the \ 2 minimization procedure. 



bution is centred on (Ci) « 120, the value obtained from 
the ROSAT XRB spectral fits in this band (eg Georgan- 
topoulos et al. 1996). The results cluster around some 
preferred values as can be seen in figure 1, in which 
we show as hatched histograms the distribution of the 
input parameters and as thick lines the best-fit output 
parameter distribution. The most probable values of the 
fitted parameters as well as their standard deviation 
over the 1000 minimizations are presented in table 1. 
For |6| > 20° we have typically that \ 2 « 33600 for 
33064 degrees of freedom and hence this model can- 
not be rejected. The Galaxy contributes a significant 
fraction (E ~ 20 — 30 per cent) of the average total 
ROSAT 1.5 keV X-ray emission. As discussed earlier 
this Galactic component could arise as contamination 
from lower energies (eg from the 0.17 keV Raymond- 
Smith Galactic component) due to the coarse energy 
resolution of the ROSAT PSPC. Although this percent- 
age is rather high it is not inconsistent with XRB spec- 
tral fits in deep ROSAT and ASCA pointings: the excel- 
lent spectral resolution ASCA spectrum of the XRB in 
the 0.4-10 keV band (Gendreau et al. 1995) presents a 
steep upturn at 1 keV, suggesting a high contamination 
of the 1.5 keV ROSAT band from lower energy pho- 
tons. The scale height of the finite emitting disk is ~ 16 
kpc, higher than both the value obtained by Pietz et al. 
(1998) in the 0.75 keV ROSAT band and by Iwan et al 
(1982) in the 2-60 keV HEAO-1 band. Using the above 
best-fit model, the predicted Galactic halo luminosity is 
L x ps 2 x 10 39 ergs The best-fit extragalactic com- 
ponent arises to ~ 100 ± 11 X 10~ 6 cts s _1 arcmin -2 . 
This value is in excellent agreement with the extrap- 
olation of the HEAO-1 3-60 keV data of Marshall et 



Table 1. Fitted Galaxy model parameters for \b\ > bn m . 

him C b Eh (kpc) R d (kpc) 

20° 95±13 0.29 ±0.06 15 ±5 27±5 

30° 100±11 0.20 ±0.05 17 ±5 27±5 

20° 120 (frozen) 0.28 ±0.10 4 ±2 15±4 

al. (1980) in the ROSAT band but ~20 per cent lower 
(~ 1.5a) than the normalization of the extragalactic 
power-law component obtained from XRB spectral fits 
in deep, high galactic latitude ROSAT fields (Hasinger 
1996). If instead we force the extragalactic contribu- 
tion to the value obtained from the ROSAT XRB fits 
(~ 120 x 10 -6 cts s -1 arcmin -2 in this band, eg Geor- 
gantopoulos et al. 1996) then the best-fit parameters are 
h = 4±2 kpc, R d = 15±4 kpc and E = 0.28±0.10 con- 
sistent with the values of Pietz et al. (1998). We finally 
note that, had we not excluded the North Polar Spur 
region in our fit, with all four parameters (Cb, E,h, Rd) 
free, we would have erroneously obtained a significantly 
higher Galactic fraction (E ~ 40 per cent) as well as a 
larger scale height (~ 19 kpc) but with an unacceptable 
fit in this case (Xred ~ 1-2). 



3 XRB ROSAT DIPOLE 
3.1 Dipole fitting procedure 

The multipole components of the ROSAT X-ray inten- 
sity are calculated by summing moments. The dipole 
moment is estimated by weighing the unit directional 
vector pointing to each 40 2 arcmin 2 ROSAT cell with 
the X-ray intensity d of that cell. We normalize the 
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dipole by the monopole term (the mean X-ray intensity 
over the sky): 

'D| _ Y^m 



V 



(2) 



M £C* 

We attempt to estimate the cosmological XRB dipole 
by applying the above procedure to the ROSAT counts 
after subtracting our best Galaxy model (table 1) and 
the regions of Galactic and "local" extragalactic X-ray 
emission. To this end we mask: (a) The Galactic plane, 

(b) the area dominated by the Galactic bulge and the 
North Polar Spur (see definitions in section 2.1) and 

(c) a region of 10° radius around the Large Maggelanic 
Clouds (we have verified that small variations in the 
limits of all the above regions do not change appreciably 
our results). 

We use two methods to model these regions: the 
first consists in substituting the observed intensity with 
the mean value estimated at high galactic latitudes (ho- 
mogeneous filling procedure) and the second based on 
a spherical harmonic extrapolation procedure (cf. Yahil 
et al f986). Since the later is slightly more involved we 
briefly review the method which is based on expand- 
ing the sky surface density field a($,tp) in real spherical 
harmonics: 



(3) 



where $ = 90° — Gal. latitude, tp = Gal. longitude (do not 
confuse the multipole I with the Galactic Longitude). In 
this formulation the normalized dipole is defined as 

1/2 



V 



3a° 



i 

E 

m — — 1 



I 1 \2 

(a m ) 



(4) 



where the factor 3 enters for consistency with the def- 
inition of eq|| The observed, £(#, tp), and intrinsic 
surface density field o~($,tp), are related according to: 
£(#, <p) = M(-d,tp) a(-d,p), where the mask M{"d,tp) 
takes values of I or depending on whether the ($, tp) 
direction points in an observed or excluded part of the 
sky, respectively. Since we are interested in recovering 
the dipole (I = 1) components of er($, tp), the correction 
terms should at least involve the quadrupole (I = 2) 
components. Expanding £($, tp) up to the quadrupole 
order and allowing for the orthogonality relation of the 
Legendre polynomials, we can express the observed co- 
efficients A™, in terms of the intrinsic ones, o™, forming 
a 9 x 9 matrix, the inversion of which then gives a™. A 
more accurate procedure would entail an expansion to 
higher order I's (cf. Lahav et al. 1994) but to recover 
a smooth underlying dipole structure (in which higher 
order Vs are negligible) we have verified, using mock 
samples, that the above procedure recovers extremely 



accurately the direction and amplitude of the true un- 
derlying dipole. 

3.2 Dipole Results 

In table 2 we present our main results for different treat- 
ment of the data. It is evident that when using the 
raw ROSAT data, the dipole points roughly towards the 
Galactic centre (in agreement with Kneissl et al 1997). 
However, when we exclude both the Galaxy and the 
North Polar Spur, the measured dipole is in much bet- 
ter directional agreement with the CMB dipole. For the 
homogeneous filling method we find <5# cm (, < 20° exclud- 
ing the Galactic plane below \b\ = 20° or \b\ = 30°. For 
the spherical harmonic method the misalignment angle 
is larger, S6 cm i, > 39°. The Virgo cluster, being quite 
near and very bright in X-rays, could contribute signif- 
icantly to the measured dipole. Excluding an area of 4° 
radius around the Virgo cluster (I, b ~ 287°, 75°), which 
is very apparent in the X-ray map, we find that it is 
responsible for ~ 20 per cent of the total dipole, while 
the dipole misalignment angle with the CMB increases 
by- 15°. 

However, it should be expected that many 
Galactic sources, probably dominating the higher 
ROSAT counts, are still present in the data and could 
affect the behaviour of the extragalactic XRB dipole. 
We therefore present in figure 2 the misalignment angle 
between the ROSAT and CMB dipoles as well as the 
normalized dipole amplitude, 23, as a function of the 
ROSAT upper count limit (C up ). The errorbars have 
been estimated by using the different Galactic model 
parameters resulting from the \ 2 fits of section 2.1 and 
presented in figure 1. We do find that our main results 
are very robust in such variations of the Galactic model. 
The two methods, used to mask the excluded regions, 
give consistent dipole results for C up ^ 140 (which cover 
~ 97 per cent of the unmasked sky). For this limit the 
ROSAT -CMB dipole misalignment angle is < 26° and 
33° for the homogeneous filling and spherical harmonics 
methods respectively. It is evident that the ROSAT - 
CMB dipole misalignment angle increases substantially 
when we include the few higher intensity cells, before 
however the Virgo cluster (which enters at C up ^ 200 
counts) starts reducing again the misalignment angle, 
as can be clearly seen in figure 2. The interpretation 
that the high intensity (C^ 140) cells are associated 
with Galactic sources is supported by the fact that when 
we include these few cells the resulting dipole direction 
moves towards the Galactic centre. We therefore con- 
sider as our best estimate of the XRB dipole its value 
at C up « 140, for which both methods used to model the 
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Table 2. The ROSAT Dipole Results for |6| > b H 





mask model 




1° 


i>° 


^CMBq 


CMB l q 


V 


sky masked 


raw counts 


Homogeneous 


20 


342.7 


20.4 


67.6 


59.5 


0.050 


36% 




Sph.Harm. 


20 


343.0 


11.0 


74.5 


63.8 


0.092 


36% 


Galaxy, North Polar Spur 


Homogeneous 


20 


280.3 


29.1 


21.2 


3 


0.011 


47% 


& Magellanic Clouds excluded 


Sph.Harm. 


20 


318.3 


20.1 


51.1 


38.5 


0.037 


47% 




Homogeneous 


30 


290.9 


41.6 


19.8 


16.1 


0.010 


63% 




Sph.Harm. 


30 


313.5 


14.9 


51.7 


36.7 


0.049 


63% 



masked areas agree and the XRB-CMB dipole misalign- 
ment angle is minimum. 

To take into account all possible sources of uncer- 
tainty, we use a Monte-Carlo simulation approach in 
which we vary all the model parameters within their 
range of validity. Using 6000 dipole realizations to take 
into account (a) the uncertainties of the Galactic model 
subtracted from the raw counts, (b) the different meth- 
ods used to mask the excluded sky regions (c) the dif- 
ferent galactic latitude limits and (d) variations of the 
excision radii around the bulge and the North Polar 
Spur we conclude that the XBG ROSAT dipole has: 



2X 



i 0.017 ± 0.008 (I, b) w (286°, 22°) ± 19° 



which deviates from the CMB dipole directions in the 



heliocentric and Local Group frames by SOcmb^ 



30° 



and <5#cmb lg ~ 10° i respectively. It is interesting that 
the ROSAT dipole is nearer to the Local Group frame 
CMB dipole direction. Our results are consistent with 
the HEAO-1 (2-10 keV) dipole (Shafer & Fabian 1983) 
which points in a similar direction (282°, 30°), albeit 
with a larger uncertainty, but has a lower amplitude: 
©heao-i ~ 0.005. 



0.04 




100 



1000 



3.3 Interpretation 

The motion of the Sun with respect to an isotropic ra- 
diation background produces a dipole in the radiation 
intensity according to: 
8C 



-j-j = (3 + a)V cos6»/c 



(5) 



where a is the spectral index of the radiation (C cx v a . 
For the 1.5 keV ROSAT band we have a ~ 0.4 (Gen- 
dreau et al. 1995). If the ROSAT dipole was totally 
due to the motion of the Sun with respect to the XRB 
(Compton-Getting effect) then we would obtain a solar 
velocity with respect to the XRB of Vq = 1300 ± 600 
km/sec, which should be compared with Vq = 369 
km/sec with respect to the CMB (Lineweaver et al 
1996). 

We therefore verify that the observed XRB 
ROSAT dipole cannot be only due to the Compton- 
Getting effect but it is significantly contaminated by the 



Figure 2. Dipole results as a function of upper count limit, 
Cup- The filled and open points correspond to the homo- 
geneous fill and Spherical harmonics mask methods respec- 
tively. The scatter correspond to the variations of the Galaxy 
model, as given by all the different x 2 fits (see section 2.1). 
The corresponding solutions, excluding the Virgo cluster, are 
shown as the continuous and dashed lines. 



dipole produced by X-ray emitting sources that trace 
the large-scale structure. We estimate the large-scale 
dipole component of the XRB by subtracting from the 
R OS A T map the expected Compton-Getting dipole and 
we obtain a dipole with £>lss ~ 0.9X>xrb pointing to- 
wards (1, 6) ~ (284°, 18°). This is in general agreement 
with Lahav et al (1997) who found that the two effects, 
contributing to the XRB dipole, are of the same order 
of magnitude. 
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4 CONCLUSIONS 

We have estimated the dipole of the diffuse 1.5 keV X- 
ray background using the ROSAT all-sky survey maps 
of Snowden et al (1995). We have first subtracted from 
the ROSAT counts the diffuse Galactic emission (the 
halo and bulge components) as well as local extended 
features such as the North Polar Spur. The Galactic 
halo model used is that of a finite radius disk model 
with an exponential scale height (eg Iwan et al 1982). 
The mean Galactic X-ray component is ~ 20 — 30 per 
cent of the background and the scale height and radius 
are ~ 16 ± 5 and 27±5 kpc respectively. We model the 
excluded regions by either homogeneously "painting" 
these regions with the mean X-ray count (derived after 
subtracting the Galaxy) or using a spherical harmonic 
expansion of the X-ray surface intensity field. 

We estimate that the ROSAT XRB dipole is point- 
ing towards (l,b) « (286°, 22°) ± 19, within ~ 30° of 
the CMB direction and having a normalized amplitude 
of ~ 1.7 per cent. The dipole direction is in agree- 
ment with previous estimates in hard X-rays (Shafer 
& Fabian 1982) but the positional errors have now been 
improved. We also find that the two effects expected 
to contribute to the XRB, ie., the Compton-Getting ef- 
fect and the anisotropy due to X-ray sources tracing 
the large-scale structure, are of the same order of mag- 
nitude, in general agreement with the predictions of La- 
hav et al (1997). However, the latter dominates having 
2?lss ~ 0.91?xrb. Finally we estimate that the nearest 
cluster, Virgo, contributes ~ 20 per cent to the total 
measured XRB dipole. 
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